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NONLINEAR  EVOLUTION  OF  LONGITUDINAL  PLASMA  WAVES 

Young-ping  Pao 
Courant  Institute  of  Mathematical  Sciences 
New  York  University 
New  York,  New  York   10012 


ABSTRACT 
Long-time  nonlinear  behavior  of  single-mode  longitudinal  plasma 
waves  is  studied  on  the  basis  of  Vlasov  equation  with  Fokker-Planck 
collision  terms.  Resonant  layer,  trapped  island,  collisional  sublayers 
and  x-point  neighborhoods  are  analyzed.  Nonlinear  evolution  equation 
is  obtained  and  wave-generated  current  is  also  evaluated. 
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I.   INTRODUCTION 

This  paper  presents  an  analytic  treatment  of  nonlinear 
longitudinal  plasma  waves.  In  the  linear  approximation,  these  waves 
decay  according  to  the  renowned  Landau  formulus  *.  Numerous  efforts 
have  been  made  to  examine  their  nonlinear  effects.  For  multi-mode 
incoherent  waves,  a  collisionless  turbulent  theory  had  been  developed 
on  the  basis  of  selective  mode  coupling,  giving  rise  to  quasilinear 
diffusion  .  At  the  same  time,  a  nonlinear  collisionless 
calculation''  for  single-mode  waves  yields  a  decay  rate  slower  than 
the  Landau  result. 

On  a  longer  time  scale,  collisional  effects  play  an  important  role 

in  the  nonlinear  evolution  of  single-mode  waves.   The  nonlinear 

q 
collisional  theory  of  Zakharov  and  Karpman7  showed  that  the  wave 

1/2 

amplitude  A(t)   decays  according  to  A  ■  -vB/A   ,   where  v  is  the 

electron  collision  frequency  and  6  a  positive  constant. 

However  there  appear  several  logical  gaps  in  the  analyses  of  Ref. 
9.  The  most  serious  gap  is  the  neglect  of  the  inertial  term  in  the 
non-resonant  region,  which  leads  to  a  direct  violation  of  the  Poisson's 
equation.  A  detailed  discussion  of  this  point  is  given  in  Section  VII. 
The  integration  constant  for  the  lowest  order  resonant  layer  solution 
in  Ref.  9  is  also  in  error,  which  directly  affects  the  final  result 
for  the  nonlinear  decay  rate. 

These  considerations  motivated  us  to  re-examine  the 
one-dimensional  problem  of  a  single-mode  electrostatic  wave  train  of 
small  amplitude  propagating  in  a  Maxwellian  background  plasma  with  low 
collision  frequency  v.  We  shall  adopt  the  more  realistic  model  of  full 
Fokker-Planck  collision  terms  and  carry  out  a  systematic  asymptotic 


-3- 
analysis  by  treating  the  resonant  and  non-resonant  regions  separately, 
Our  result  leads  to  the  nonlinear  equation 


A(t)  ■=  -v(YA  ♦  aA1/2)  (1) 


where  the  linear  decay  constant  1  is  positive  and  the  constant  a  is 
positive  for  z  <  12  and  negative  for  z  ^  13.  z  being  the  charge  number 
of  the  ions. 

The  linear  term  YA  represents  the  damping  due  to  collisions 
between  the  bulk  thermal  electrons  and  the  ions  while  the  nonlinear 
a-term  arises  from  the  collisional  re-distribution  of  the  space  charge 
of  the  resonant  electrons. 

The  main  difference  between  the  present  result  in  eq.   (1)  and 

I/O 

that  of  Ref.   9  is  the  A    dependence  of  the  nonlinear  term  in  (1) 

-1  /2 
while  Ref.   9  gives  a  nonlinear  term  with  A     dependence.  The 

coefficient  a  in  (1)  is  also  completely  different  from  that  given  in 

Ref.   9. 

The  basis  of  th  present  analysis  is  the  Boltzmann-Vlasov  equation 

with  Fokker-Planck  collision  terms.   Nonlinear  and/or  collisional 

effects  become  important  in  the  velocity  space  resonant  layer  where  the 

electron  longitudinal  velocity  approaches  the  wave  phase  velocity.  The 

dynamics  of  the  resonant  electrons  can  be  dominated  by  nonlinearity  or 

collisions  depending  on  the  magnitude  of  the  parameter  A  *=  v/tJ  where 

2 
e   is  the  magnitude  of  the  wave  amplitude.   For  large  A,  collisions 

prevail  over  nonlinearity  and  therefore  linear  theory  applies  and  leads 

to  Landau's  result.   On  the  other  hand  when  A  is  small,  nonlinearity 
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overtakes  collisional  effects  for  the  resonant  electrons  and  the  result 
Is  qualitatively  different  from  the  linear  theory. 

In  this  paper,  only  the  more  interesting  nonlinear  case  of  small  X 
is  considered.  In  this  case,  although  nonlinearity  dominates  the 
resonant  electron  dynamics,  collisional  effects  also  play  an 
indispensible  role  in  the  long  time  behavior  of  the  wave  train.  The 
importance  of  collisions  comes  about  in  two  ways.  First,  acting  as 
higher  order  corrections,  the  collision  terms  determine  the  lowest 
order  distribution  function  of  the  resonant  electrons.  Secondly, 
collisions  provide  the  necessary  transition  between  the  trapped  and 
circulating  electrons. 

The  main  nonlinear  feature  of  the  problem  is  that  the  deviation  of 
the  resonant  electron  distribution  from  the  Maxwellian  is  of  order  e 

which  is  proportional  to  the  square  root  of  the  wave  amplitude,  while 

2 

in  the  linear  theory  the  deviation  is  of  e  order. 

The  basic  equations  and  variables  are  given  in  the  next  section 
and  solution  to  these  equations  for  the  non-resonant  electrons  are 
obtained  in  Section  III. 
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II.   EQUATIONS  AND  WAVE  VARIABLES 

Consider  a  single-mode  wave  train  with  wave  number  tc  and  basic 

frequency  u>.   By  normalizing  the  distance  x  with  1/k,  the  time  t  with 

1/2 
1/tcvef  where  vg  ■  (Te/m)    and  te  is  the  electron  temperature  at 

t  -  0.  We  introduce  the  wave  variables  6  and  the  phase  shift  4>(t)  as 

6  =  x  -  ct  +  4»(t)  , 


(2) 


C  «=  u/KVg  , 


and  express  the  electric  field  E  as 


c  .  .  <mve  dp 

e   de  ' 


where  e  is  the  (signed)  electron  charge.  We  shall  also  normalize  the 

electron  distribution  function  F  by  NZ/v|,  N  being  the  ion  number 
density. 

In  terms  of  these  variables,  the  basic  equations  become 


—  *    (v  -c)  — -  -  -I  - —  =  vC(F,F)  (3) 

dt  x    96   d6  9v 


£♦-  (-^-)2(l  -  JFd3v)  ,  (4) 

de2    KVe      J 


where  the  ions  are  assumed  to  have  a  constant  Maxwellian  distribution, 
and  the  electron  plasma  frequency  we  and  the  electron  collision 
frequency  v  are  given  by 
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u>|  -  4irNZe2/m  , 


v  -  2TTNZC14  In  A/<m2Vg  , 


A  Is  the  Coulomb  logarithm  and  the  detailed  form  of  the  Fokker-Planck 
collision  terms  C(F,F)  can  be  found  in  Ref.   10. 

Because  of  azimuthal  symmetry,  the  electron  distribution  function 
F  has  no  azimuthal  angle  dependence  and  may  be  written  as 

F  ■=  F(v  v  ,6,t)  , 
1 

pop 

where  x  is  the  direction  of  propagation  and  v ,  »  vr;  +  v^. 

1    Y 
Since  the  wave  amplitude  is  small,  the  potential  $  may  be  expanded 

as 


4  •=  e2<t>0(e,t)  +  e\(6,t)  +  ...  (5) 


where  the  lowest  order  form  is  <J>0  =  A(t)  cos  6  as  anticipated  from  the 
linear  theory  with  the  provisor  that  A(t)  varies  with  t. 

The  mathematical  problem  is  to  solve  (3)  and  (k)  for  c  <<  1  and 
v  <<  1  .  For  sake  of  simplicity,  it  is  also  assumed  that  c  >>  1  in 
conformity  with  many  wave  experiments  where  the  wave  phase  velocity  is 
much  larger  than  the  electron  thermal  velocity  vg. 

To  solve  the  Boltzmann  equation  (3)i  we  divide  the  velocity  space 
into  two  regions:  a  non-resonant  region  with  |vx~cl  >>  ei  where  the 
linear  theory  provides  the  lowest  order  solution  and  a  resonant  layer 
with  vx  -  c  -  ei  where  the  electrons  must  be  treated  nonlinearly.  The 
scale  parameter  e1  will  be  made  explicit  in  the  next  section.   Once  the 
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solutions  are  obtained  for  the  two  regions,   they  can  be  substituted 
into  the  Poisson's  equation   (M)   and  the  result  may  be  symbolically 
expressed  as 


2 

14  '   -k2*  +  Pi(e,A.A.*.c)  ,  (6) 

de^ 


7 
where  the  K  -term  is  the  linear  theory  result  from  the  lowest  order 

non-resonant  solution  and  p<  is  the  higher  order  contributions  from  the 

non-resonant  as  well  as  resonant  electrons. 

The  lowest  order  result  of  (6)  is 


d2*o     2 

de^ 


which  leads  to  the  linear  dispersion  relation 


K2(k,c)  ■=  1 


and  our  lowest  order  result   for   $: 


<}>0  =   A(t)    cos   6 


In  the  next  order,  we  have 


d2^ 


♦l  *  Pi 
de^ 
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It  then  follows  that  p.  must  satisfy  two  compatibility  conditions, 
namely 


\2* 


<J6  p1  (e,A,A,\J/,c)  sin  6=0, 


(7) 


2v 


de  p1  (e,A,A,4>,c)  cos  e  =  o 


(8) 


and  these  two  equations  determine  A  and  ij>  and  therewith  the  evolution 
of  the  wave  train. 

On  the  other  hand,  the  solutions  for  F  in  the  non-resonant  and 
resonant  regions  can  also  be  used  in  the  energy  equation  to  evaluate 
the  change  of  background  Kaxwellian  in  time,  which  will  be  done  in 
Section  VII. 

The  details  of  the  last  two  equations  will  be  postponed  until 
Section  VI;  for  the  moment  we  shall  turn  to  the  non-resonant  region 
analysis . 
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III.   NON-RESONANT  REGION 

Solutions  for  the  non-resonant  region  can  be  obtained  by  iteration 
from  the  linear  theory  results.  The  distribution  function  F  can  be 
decomposed  as 


F  =  Fn^yV.'t)  +  Hv  .v.e.t)  , 
1  1 


(9) 


where  F   is  the   6-averaged  part  and  f  the  6-dependent  part  with  the 
6-average 


<f>fl  r  ^~ 
6  -  2tt  J 


r2n 


f(vv,v, ,e,t)de  =  o  , 
1 


do) 


and  f  is  periodic  in  e  with  period  2ir .  The  t-dependence  in  FQ     and  f 
allows  for  slow  variation  in  time. 

Two  separate  equations  can  be  derived  for  FQ  and  f.  First, 
equation  (3)  can  be  e-averaged  with  the  result 


—  -   <£*  —> 
3t     d6  3vx  e 


vC(F0,F0)  +  v<C(f,f)> 


(11) 


which  describes  the  time  evolution  of  F.  Then  subtracting  (11)  from 
(3)  yields  the  following  equation 


(v  -c)  |£  -  *! !!£  -  vc,(f)  -  <?±1L>6-*L 

x  36        d6   3vx  L  de   3vx   e        3t 


(12) 


for  f,  where  C,  is  the  linearized  collision  terms 
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CL(f)  -  Cee(Fo'f)  +  <Wf'Po>  +  Cei(f)  • 

In  the  non-resonant  region  | vx-c |  >>  e1 .  The  last  two  terms  in 
(12)  will  be  shown  to  be  of  v  order.  Therefore,  if  z*  >>  v,  then  all 
the  right  hand  side  terms  of  (12)  are  much  smaller  than  its  left  hand 
side  terms.  This  then  calls  for  the  expansion 


f  "  fo  +  f1  +  ' 


?        ? 
where  f  -  e  ,  f.,  -  ve  and  they  satisfy,  respectively, 


,     .  3f0    6, 9<w    . 


(13) 


3f ,  3fn  rtA    3fn 


(11) 


Note  that  in  eq.  (13)  we  have  carried  the  higher  order  term 
3f0/9vx,  just  to  save  iteration  steps.  In  fact,  integrating  (13)  along 
the  characteristic 


Y  =  I  (vx-c)2  *  $ 


yields  the  result 


d$ 


3F, 


de  (2Y-2$)1/2  8vx 


de 


(15) 
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and  for  $  small,  the  lowest  order  f  is  just  the  linear  solution 


3F 
f  -  f   -  _^_   °  .  (16) 

o    oo   vx  -  c  8vx 


It  is  clear  from  (15)  that  fQ  -  f0^vx'v  '^'l^  t0  a11  orders  of  <t> 
and  consequently  the  last  term  in  (1M)  vanishes.  Then  the  lowest  order 
solution  of  {^k)   becomes 

f1  =  e2(f1s  sin  6   f1c  cos  6)/(vx-c)  ,  (17) 

fls  =  A(t)vCL(MF0)  -  A(t)MF0  ,  (18) 

f1c  =  -A(t)*(t)MF0  ,  (19) 

MF0  =  OF0/9vx)/(vx-c)  .  (20) 

Now  we  express  F  in  the  non-resonant  region  as 


F0  ■=  FM  *  F1  +  ...  ,  (21) 


where  FM  is  the  Maxwell ian  which,  in  our  normalized  variables,  becomes 
FM  =  (2^T)"3/2  exp  (-v2/2T)  , 

where  T  -  T(t)  is  the  electron  temperature  at  t  normalized  by  the 
electron  temperature  Te  at  t  -  0. 

The  deviation  F1   from  the  Maxwellian  is  responsible  for  the 
generation  of  current  drive  and  is  therefore  a  quantity  of  interest. 
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Now  substituting  (21)  into  (11)  yields  the  equation  in  the  non-resonant 
region: 


vCL^i)-^-<^^>e-v<c(f0.f0)>  , 


(22) 


which  can  be  solved  for  F1  if  3FM/3t  is  known. 

To  evaluate  the  rate  of  change  of  FM,  we  return  to  the  whole 

velocity  space  and  use  the  original  equation  (3).   Multiplying  (3)  with 

2 

mv  /2  and  integrating  over  6  and  v  yield  the  energy  equation 


Imf 
2 


■-I 


d*     ■}         i 
mv   <— I  F>c,6Jv   *   -  v 
x   de    6   -   2 


mv^Cei(FM)dv  , 


where  the  last  term  represents  the  energy  transfer  between  electrons 
and  ions  and  is  of  the  order  v(m/mi),  n^  being  the  ion  mass.  The 
ion-electron  energy  transfer  can  be  evaluated  explicitly  and  the  final 
form  of  the  energy  equation  then  becomes 


It 

2 


u  =  ^z(m/mi)/(2iTT3)1/2  , 


(23) 


where  T^  is  the  ion  temperature  normalized  by  Tg. 

It  must  be  emphasized  that  the  integral  in  (23)  covers  both  the 
resonant  and  the  non-resonant  regions.  We  now  proceed  to  the  resonant 
layer  analysis  where  nonlinearity  plays  an  important  role. 
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IV.   RESONANT  LAYER 

In  the  resonant  layer,  we  introduce  the  new  variable  u  and  write 


vx  -  c  -  Clu 


where  e1  is  the  length  scale  in  the  layer.  The  Vlasov  equation  (3)  can 
then  be  written  as 


3F   e2A  sin  6  3F    .  ._,   3F  ._.,. 

e.u  — -  +  —  -  vC,  (F)  -  _  ,  (2U) 

1   36      t1     3u     L      3t 


and  the  linearized  collision  terms  C«  (F)  takes  the  simplified  form 


CL(F)  =  (I*,)  *  [ilZgfl  |I]  ♦  4  »  (1  |£  ♦  F)  (25) 

b  3y    v3   3p    v2  3v  v  3v 


in  terms  of  the  pitch  angle  variable  y  ■  vx/v* 

For  resonant  electrons,  most  of  the  collisions  are  with  the  bulk 
thermal  electrons  and  ions.  Since  the  longitudinal  velocity  of  the 
resonant  electron  is  close  to  c  which  is  much  larger  than  the  thermal 
speed  of  the  bulk  electrons  and  ions,  the  Fokker-Planck  collision  terras 
can  be  simplified  to  yield  the  result  in  (25).  In  terms  of  the 
variables  u  and  v  ,  eq.   (25)  can  be  further  simplified  as 


CL<F>  -^CR(F) 
cfc^ 


where  the  resonant-layer  collision  terms  CR(F)  are 
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CR(F).^4+E1CL1^+^V  . 

v  3u2        9u 

v0  -  v[(1+z)v2  ♦  2T]  , 

1  (26) 

L1  -  -2[z  +  (1  +  z)D  ]  , 

L?  -  (1+z)v~2D2  ,      D,  =  v  —   , 
2  11        1    1  3v 


and  v  is  considered  to  be  of  order  unity. 

Using  e^c  as  a  small  parameter,  we  obtain  the  lowest  order  result 
of  (24)  as 


,,  3F  w  c2A  sin  e  3F    vo  32F  ,__. 

36      e2    9u   c3c3  3u2 


where  we  have  neglected   the   inertial  term  in  (24)   by  assuming 

2  o 
3/3t  <<  v/e^c-*;   i.e.,   we  look  for  long-time  behavior  of  the  wave 

propagation  and  the  neglect  of  the  inertial  term  will  be  justified 

2  3 

a  posteriori.   In  the  opposite  limit  of  3/3t  >>  v/ztc-',   the  collision 

terms  can  be  neglected  in  (24)  and  the  problem  reduced  to  the 

7 
collisionless  nonlinear  case  that  has  already  been  treated  by  O'Neil. 

The  nature  of  eq.   (27)  depends  on  the  parameter  X  -  v/e-^c  .   For 

large  X,  collisions  dominate  nonlinearity  while  small   X  means  that 

nonlinearity  prevails  over  collisions.   We  discuss  these  two  cases 

separately. 

(i)  X  >>1.   We  choose  z,    ■  v  ^c   and  rewrite  (27)  as 
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[0*z)v2   ♦  2T]  tl  -  u  E  -    X_2/3A  sin   6  II   , 
1  3u2  3e  3u 


—2/  3 

which  can  be  solved  by  expanding  in  the  small  parameter  X   J  and  the 

lowest  order  solution  can  be  constructed  from  Airy  functions.  The 
final  result,  obtained  by  matching  the  resonant  layer  solutions  to  the 
non-resonant  solutions,  agrees  with  the  linear  theory  of  Landau 
damping.  Although  collisions  dominate  over  nonlinearity ,  we  must  bear 
in  mind  that  v  <<  1  ;  i.e.,  collision  frequency  is  small  in  comparison 
with  kc.  It  is  not  surprising  that  this  should  lead  to  the  result  of 
collisionless  linear  theory. 

(ii)  X  <<  1  .   In  this  case,  we  choose  z-^    «  e  and  rewrite  (27)  as 


uE   ♦  A  sin  6  E   -  A.lfl  ,  (28) 

36  Su    o  Su2 


where  we  have  XQ  -  AvQ/v  <<  1 . 

Our  main  interest  lies  in  Case  (ii)  where  nonlinearity  plays  a 
dominant  role.  For  our  purpose,  the  higher  order  collision  terms  in 
(26)  are  also  important  and  we  therefore  shall  retain  these  terms  in 
(28)  and  consider 


u  E  *   A  sin  6  —  -  XCR(F)  ,  (29) 

96  3u     R 


where  Cn(F)  contains  all  the  terms  in  (26) 
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We  can  now  solve  (29)  by  expanding  for  small  A: 


F  "  Fo  *  XF1  + 


and  the  lowest  order  equation  for  F  is 


3fn 


where  the  characteristic  variable  £  is  just  the  one-dimensional  energy 
integral 


and  the  solution  is 


5  -   u2  +  A  cos  e   ,  (30) 


Fo  =  Fo^> 


The  constant  %     contours  in  the  (u,6)-plane  form  an  island 

structure  where  the  separatrix  at  £  =  A  separates  the  circulating 

electrons  for  £  >  A  from  the  trapped  electrons  for  £  <  A;  see  Fig.   1 . 

The  solutions  take  different  form  for  the  two  classes  of  electrons  as 
shown  in  the  following. 
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IV. A.   Solutions  for  Circulating  Electrons 

The  first  order  terms  In  (29)  satisfy  the  equation 


SF^.e)   x 


de 


-  h^  • 


(3D 


and  since  F.   must  be  periodic  in  6,  we  can  integrate  (31)  over  6  for 
fixed  £  to  obtain 


<CR(F0)/u>  =  0  , 


(32) 


where  the  bracket  denotes  the  average 


<F>  =  _ 

2ir  J 


2n 


FU,6)d6  ,      £  >  A 


(33) 


with  £  kept  constant;  this  is  different  from  the  e-average   in   (10) 


which  keeps  vx  fixed, 


We  now  solve   (32)   for  F0  by  substituting  (26)  into  (32)  and 
expanding  in  the  small  parameter  e.|C,  where  now  e1  -  c.   If  we  express 


F0  as 


F0U)  -  Fc  +  ecG(0  +  e2c2H(0  ♦ 


(3*0 


where  F  is  the  Maxwellian  distribution  at  vx  =  c,  viz., 


Fr   -    (2t»T)  3/2  exp  [-(v2+c2)/2T] 

1 


Then  substitution  of  (26)  into  (32),  up  to  c2c2  order,  yields 
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u  *>  2 
3u 


>  -  0  , 


(35) 


vo  J  *2*± 


<- 


1     98+     l 
>  +  <_  L,  -^>  +  <-  L5F_>  -  0  , 


v   u  gu2      u   >  9u 


u 


'2l  c' 


(36) 


where  we  have  used  the  +  and  -  subscripts  to  denote  the  directions  of 
the  electron  motion  with  u  >  0  and  u  <  0,  respectively. 

As  e;  ■»  <*,   or  u  +  ±» ,  FQ  must  match  the  solution   in   the 
non-resonant  region  for  v  +  c,  namely 


2  2 
_    _  r,    ec     EC   2 

i       2T^ 


]  • 


This  then  shows  that  the  boundary  conditions  for  G  and  H  are: 


G±U)  -  ♦  /2l   Fc/T  , 
H.(5)  -  CF./T2  , 


(37) 
(38) 


for   u  ■*  ±" . 

The     solution     for      (35)   subject   to  the  boundary  condition   (37)    is 
given  by 


G+U)    =    *(Fc/T)[/2~l  -    X(0]    . 


X(0 


[^rUr  "  —  ]dn  . 


5     <|u|>       /^ 

|u|    -    (2n  -   2A  cos  e)1/2 


(39) 
<uo) 
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To  solve  (36),  we  use  3/3u  -  u3/H  and  rewrite  eq.   (36)  as 


[_£  <u>H^  *  L1G±  ♦  <u>L2Fc]'  -  0  ,  (41) 


where  the  prime  denotes  differentiation  with  respect  to  £. 

We  now  integrate  (41)  and  set  the  integration  constant  to  zero  on 
account  of  (38)  and  the  final  result  becomes 


H+U)  -  FCT~2[1  +  P(OQ(v^/2T)/<|u|>]  ,  (42) 

P(0  =  /2i  -   <|u|>  -  xU)   ,  (13) 


Q(x)  =  (2x-oz)/(x+o)  ,     0  =  —1-.  ,  (44) 

(1+z) 


where   P(£)  ■»  0  as  u  ■♦  ±»  so  that  (42)  satisfies  (38)  and  x(0  is  given 
in  (40). 

IV. B.   Solutions  for  Trapped  Electrons 

For  trapped  electrons  with  I,  <  A,  the  Z,  contour  is  bounded  in  6 
between  e^  to  the  left  of  6  =  it  and  ep  to  the  right  of  e  =  tt  ,  0£  and  ep 
being  the  two  roots  of  cos  6  -  £/A. 

Similar  to  the  circulating  electron  solutions,  the  lowest  order 
(in  A)  solution  for  the  trapped  electrons  has  the  form 


Fo  -  r0U) 


and  the  A-order  equation  is  again   (31).   Now  we  cannot  use  the 
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periodicity  argument.   Instead,  we   invoke  the  continuity  of  F1   as 
u  *  ±0  and  the  result  becomes 


<C#(F0)/|u|>  -  0  , 


(45) 


where  C«  is 


C.CF)  =  ^  |u|  ±   [|u|  |I)  ♦  ccL^uj  »4  e2c2L2F  , 


(46) 


and  the  bracket  in  (45)  is  defined  as 


<F>  =  _ 

2n  J 


r  r 


F(^,e)de  , 


(47) 


for   £  <  A.    In  other  words,  for  £  >  A  the  bracket  <F>  is  defined  in 
(33)  and  for  £  <  A  it  is  defined  in  (47). 

Again,  we  shall  solve  (45)  by  expanding  in  ec: 


F0  =  Fc  -  ecG.U)  *  c2c2H»(0  * 


(48) 


which,  upon  substituting  into  (45),  leads  to  the  result 


i   _  i 


[<  u  >G«U)]   -  0  , 


(45) 


[—  <|u|>H»  4  L.G,  +  <|u|>L2F  ]' 


2'c- 


(46) 


The  solution  for  (45)  is 
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G«(0  -  c1  ♦  c2  J  dt1/<|u|>  (47) 


which  must  match  (39)  as  £  ■»  A  for  u  >  0  and  u  <  0.  This  is  impossible 
since  G  in  (39)  have  different  signs.  The  way  out  of  this  dilemma  is 
to  allow  for  a  singular  layer  near  the  separatrix  at  £  ■  A  where  the 
collision  terms  cannot  be  treated  as  small.  These  collisional 
sublayers  will  be  examined  in  the  next  section.  One  of  the  results 
from  the  collisional  sublayer  analysis  shows  that  the  constants  in  (H7) 
must  vanish,  i.e.,  c1  =  Cp  ■  0,  and  consequently 

G*U)  5  o  . 

We  can  now  integrate  (46)  with  the  result 


H*(0  =  FCT~2[1  -  Q(v*/2T)]  ,  (48) 


where  Q  and  o  is  given  in  (44).  In  deriving  (48),  we  have  used  the 
boundary  condition  that  H*  is  regular  at  the  island  center  £  =  -A. 

It  is  also  noted  that  (48)  does  not  match  (42)  at  the  separatrix 
so  that  H*  must  also  undergo  a  rapid  transition  through  the  collisional 
sublayers.  However,  H»(£)  can  match  H  (£)  at  £  =  A  by  adjusting  the 
constant  of  integration;  i.e.,  H+(£)  and  H*(£)  are  continuous  at  E,   ■  A. 

Finally,  we  note  that  we  have  neglected  the  9F/3t  term  in  (24)  in 
obtaining  the  resonant  layer  solutions.  This  term  may  be  expressed  as 


3F    ;  3F    •  3F    •  3F 

—  =  A  ♦  iJ;  —  +  T  —  . 

3t     3A     36     3t 
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The  neglect  of  the  A  and  T  terms  can  be  Justified  since,  as  we 
shall  see  later,  A  -  v  and  T  -  ve  which  is  much  smaller  than  the  XCR 
term  in  (24).  On  the  other  hand  it  will  be  shown  that  i>  ~  EC  and 
therefore  the  ^  term  can  be  neglected  in  comparison  with  the  lowest 
order  term  u9F/96.  Furthermore,  retention  of  the  \\>  term  as  a  higher 
order  term  in  (24)  or  (31)  will  not  affect  the  results  in  (32)  or  (45) 
because 


ae  u     °  de  u 


on  account  of  the  antisymmetry  in  e  and  this  justifies  the  neglect  of 
3F/3t  in  the  resonant  layer  analysis. 
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V.   COLLISIONAL  SUBLAYERS  AND  x-POINT 

We  now  proceed  to  analyze  the  collisional  sublayers  and  to 
substantiate  our  previous  claim  that  the  constants  c.|  and  C2  in  (^7) 
vanish.  Our  analysis  follows  the  method  developed  in  Ref.  11  for 
hydrodynamic  stability  problems;  for  a  somewhat  different  version  of 
this  method  see  Ref.   12. 

There  are  two  branches  of  the  sublayer  at  £  -  A,  an  upper  branch 
for  u  >  0  and  a  lower  one  for  u  <  0.  In  the  sublayer,  £  »  A  and  we 
therefore  introduce  a  new  variable  n  through 


t.    -    A  *  Aj'2n 


f0  =  fc  ♦  g0(n.e)  ♦  xj/^di.e) 


The  velocity  variable  u  can  be  expanded  as 


u  =  ±2A1/2(sin  1  e)  (1  *  \\/2x\/2h   sin2  -   e)  ,         (^9) 


provided  that  6  >>  A^7 ' U   and  2n-e  >>  \]J    .   Near  the  x-points  where 
6  -  Aq   or  2n-6  -  Aq   ,  a  different  analysis  is  needed. 

Substituting  the  expansion  for  F0  and  (il9)  into  (29)  and  using  the 
dominant  term  in  CR,  we  obtain 


3gQ   3?So  .  (50) 


'o 
36 


3n2 
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4-  (n^2)    .  (5D 


Se+        3n2    "   ,„   ..2  1    „  3n  ln    3n 


3n^        ijA  sin^  _  e 


where  the   6      variables   are 


6±   ■=    UA1/2(l    -   cos  I  e)    ,  (52) 


and  the  plus  and  minus  subscripts  denote  u  >  0  and  u  <  0,  respectively. 
As  n  •*  ±",  g0  and  g-  must  match  (39)  and  (Ml).  For  this  purpose, 
we  rewrite  g0  as 

go  =  1  G±(A)[l  +  erf(n/ze]./2]  ♦  1  cjl  -  erf(n/2e]/2]  ♦  T(n.e±)  ,   (53) 

where  the  error  function  approaches  ±1  as  n  ■»  ±°> . 

Using  a  generic  variable  E  for  6+  and  6_ ,  we  have 


9Y(n,5)  =  32Y(n,E)  (5iJ) 

9=    "    Sn2    ' 


and  Y  tends  to  zero  as  n  *  ±". 

As  e  ■*  0  and  6  *  2ir,  the  solution  in  (53)  must  match  the  x-point 
solutions  which  will  be  shown  later  to  be  a  function  of  n  only.  This 
matching  is  then  equivalent  to  matching  the  expression  in  (53)  to 
itself  at  6  -  0  and  e  -  yielding  four  conditions  on  Y±: 

For  n  >  0  , 
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Y±(n.O)  +  G±(A) 


1  G±(A)E0  +  1  clEl  ♦  Y±(n,H_)  .     (55) 


For  n  <  0 


Y±(n.O) 


1  G_(A)E0  +  1  clEl  *  Y_(n,r0)  . 


(56) 


where  G+(A)  is  given  in  (39)  and  E^    and  E2  are  given  by 


1/2, 


■n  =  En(n)  =  1  ♦  M)nerf(n/zE*/z)  , 


8A 


1/2 


)+(e  ■=  2u)  «=  e_(e  =  o) 


(57) 


The  derivation  of  (57)  and  (58)  is  similar  to  that  used  in  Ref. 
11.  Equations  (55)  and  (56)  can  be  converted  into  a  simpler  form  in 
terms  of 

S(n,=)  =  Y+(n,H)  ♦  Y_(n,E)  , 
D(n,5)  =  Y+(n,5)  -  Y_(n,5)  , 

and  the  resultant  equations  become 


S(n.O)  -  S(n.H0)  =  c^sgn  n)E1  (  |  n  | )  , 
D(n.O)  -  (sgn  n)D(n,50)  -  -G±(A)E1 ( | n  | ) 


(58) 
(59) 


Now  we  can  show  that  c1  ■  0  by  using  the  transform 


SU.5) 


S(n,=)  exp  (inc)dn 
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which  exists  and  is  continuous  since  S  *  0  as  r\  ■*■   ±».   Since  S 
satisfies  (5^),  the  transform  S  has  the  property  that 


S(c,=)  -  S(c,0)  exp  (~c2H)  ,  (60) 


while  the  transform  version  of  (58)  is 


SU,0)[1  -  exp  (-;2H0)]  -  c^U)  ,  (61) 


where  E.|    is  the  Fourier   transform  of    (  sgn  n )E1  (  | ti  |  ) . 
In  the   limit   of   ;   +   0,    eq.      (61)    becomes 


S(0,0);2Eo  -  i  c1   J    |Ti|E1(|n|>dn 


and  we  readily   conclude   from  the   continuity  of  S  that 


C.    =    0 


as  claimed  in  the  last  section.  A  direct  consequence  of  this  is  that 
SU.O)  ^  0  and  therefore 

S(n,H)  5  0  .  (62) 

Equation  (59)  can  be  solved  by  Wiener-Hopf  method,  provided  the 
singularity  of  the  half-range  transforms  at  z.  -  0  is  taken  care  of. 
The  solution,  together  with  (62)  yield  the  solutions  Y±(n,H),  but  these 
will  not  be  needed  for  our  calculation. 

Similarly,  we  can  use  the  same  procedure  for  g^  and  prove  that 
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c2  -  0 


It  is  useful  to  observe  that 

gjCn.e,  u  >  o)  =  -g^n^n-e,  u  <  0)  ,     1-0,1,     (63) 

which  follows  from  (50),  (51)  and  the  condition  that  g0  and  g1  must 
match  (39)  and  (41).  Equation  (63)  expresses  the  skew-symmetry  for  the 
sublayer  solutions. 

Finally,  we  consider  the  neighborhoods  near  the  x-points  at  u  =  0 
and  6=0,  2n.  In  these  neighborhoods,  6  -  A1M  (or  2tt-6  -  X1/lJ)  and 
u  -  X    ,  and  consequently  the  lowest  order  x-point  equation  is 


u  _—  -   A  sin  6  — -  «=  0  , 
36  8u 


which  shows  F  is  a  function  of  £  (or  n)  only,  as  claimed.  Higher  order 
solutions  can  be  constructed  but  we  shall  not  go  into  further  details 
since  they  are  tedious  and  not  needed  for  our  purpose.  Suffice  it  to 
say  that  the  sublayers  and  the  x-point  neighborhoods  provide  a  smooth 
(though  rapid)  transition  for  the  solutions  of  the  resonant  electrons. 
What  is  needed  for  our  later  calculation  are  the  conclusions 
c.|  =  c2  ■=  0  and  the  skew-symmetry  in  (63). 
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VI.   SPACE  CHARGE  AND  NONLINEAR  EVOLUTION 

To  determine  A  and  \p  from  (7)  and  (8),  it  is  necessary  to 
calculate  the  space  charge  in  the  Poisson's  equation.  We  first 
consider  the  sin  6  component  ps  of  the  space  charge  and  split  it  into  a 
non-resonant  part  pcr,  and  a  resonant  part  pQ„  as 

oil  9 1 


s 


dedvF  sin  e  =  psn  -  Psr  ,  (6*0 


P_n  -     de  sin  e    dvF  ,  (65) 

sn   Jo  -"n  " 

r2n 


Psr 


de  sin  6    dvF  ,  (66) 

J  r 


where  the  psn  integral  for  non-resonant  electrons  with  subscript  n  is 
over  |vx~c|  >  6  ancl  the  Psr  integral  for  resonant  electrons  with 
subscript  r  is  over  | v  — c |  £6,6  being  a  parameter  satisfying  6  <<  1 
and  6  >>  e. 

For  the  non-resonant  integral  in  (67),  we  observe  that  only  f,  in 
(17)  has  a  non-vanishing  contribution  while  the  integrals  over  F  and  f 
vanish  due  to  symmetry  in  6.  Substituting  (17)  and  (18)  into  (65)  and 
expanding  for  c  >>  1 ,  we  have 


psn  -   2e2c"3(A  ♦  >veA)  ♦  6-terms  ,  (67) 

Y  -  n/3(2ir)1/2  ,     ve  -  vT~3/2  ,  (68) 


where  the  6-terms  are  terms  that  depend  on  6.  Note  that  the  first  two 
terms  on  the  right  hand  side  of  (67)  respresent  the  contribution  of  the 
thermal  electrons  and  ions  while  the  6-terms  come  from  the  boundary  at 
vx  «  c  +  6  and  vx  ■=  c  -  6.   These  6-terms  must  cancel  the  boundary 


terms  in  the  resonant  part  p 


sr' 
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The  simplest  way  to  derive  (67)  is  to  substitute  (18)  into  (65), 
expand  for  large  c  and  ignore  the  gap  c-6  £  vx  <  c+6.  In  a  more 
rigorous  way  of  deriving  (67),  we  should  subtract  the  singularity  of 
the  integrand  in  (67)  as  v  ■*■  c  and  then  take  the  limit  of  large  c,  but 
this  is  very  tedious  and  we  shall  not  dwell  on  it  here. 

Now  to  calculate  (66),  we  rewrite  psp  as 


Psr  '  2* 


v  dv  o 
o  1  |Psr  ' 


'sr 


r2i 

c     d6  sin  e 
Jo 


-A 


du  F  , 


(69) 
(70) 


where  6  ■  6/t  >>  1  ,  and  we  have  used  vx~c  ■=  eu. 

We  now  change  the  variables  (u,6)  to  (£,6)  where  £  is  defined  in 

(30)  and  use  a  to  denote  the  sign  of  u;  i.e.,  a  ■  ±1  for  u  >  0  and 
u  <  0,  respectively.   Then  (70)  becomes 


psr  =  E 


r   f2lT  fA« 

I    \        de  sin  6     d£  F(S,6;a)/U 


1   ? 


e»  -  a  cos  e  , 


U  -  u  , 


(71) 
(72) 


where  the  argument  a  in  F  indicates  that  F  is  different  for  u  >  0 
(a  -  1)  and  u  <  0  (a  ■=  -1),  and  the  summation  in  (71)  is  over  a  -  1  and 
a  -  -1  . 

Noting  that  A  sin  e/U  ■=  911/96,  we  can  integrate  (71)  by  parts  in  6 
and  write  the  result  as 


psr  -  -eA~1  I  J  de  J  d£  U 


3F 

96 


(73) 
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The  e-derivative  in  (73)  is  just  the  left  hand  side  term  of  (29) 
except  for  the  different  signs  of  U  and  u.  Consequently  using 
a  -   sgn   (u),    we  may  write 


'sr 


:A   1    I   a    I 


r 


d£   CR(F) 


(74) 


and  if  we  substitute  CR  in  (29)  into  (7*0  and  use  the  E;  variable,  eq. 
(7^)  becomes 


sr 


f   de   p    [X    »    („*   »  -  F)   ♦   ecL,U  »  -   c,cVL2F]    •      (75) 
Jo  Jf*        v     3E  3E  '35  c 


AeA  '    I    \        de 

'o     JE»    v   3E 


The  first  term  in  the  bracket  of  (75)  is  a  total  derivative  and 

can  be  immediately  integrated.  At  the  upper  limit,  F  approaches  the 

Maxwellian  and  therefore  the  contribution  vanishes  on  account  of  the 

sign   in  a   except  for  higher  order  6-terms.   At  the  lower  limit  E  «=  E* , 

2 

u  «  0  and  therefore  the  u  -term  vanishes  while  the  F  terms  cancel  out 

since  I   aF(E,a)  =  F(E*,1)  -  F(E„,-1)  -  0 


as  F  must  be  continuous  at  u  ■  0. 

We  now  substitute  (75)  into  (74)  for  p 
upon  integratig  over  v  ,  and  consequently 


The  L0F  term  vanishes 
sr       d 


psr    ■=    -2ttAe2cA    1 


0  1  1 


2n     (A#      - t 

de     dEL,UF^  6-terms, 


(76) 


where  the  solution  F_  has  the  form  in  (34).  The  F  term  in  (34)  makes 
no  contribution  to  the  integral  in  (76)  since  it  is  independent  of  E. 
contribution  of  the  G  term  also  vanishes  since  G+  is  anti-symmetric  in 
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the  resonant  layer  outside  the  collisional  sublayer,  while  inside  the 
sublayer  the  solution  is  skew-symmetric;  see  (63).  The  final  result  is 
then 


Psr  -  -21tXeJ4c3A  1  I 


v  dv  L 

o  1  1 


2tt 


de 


d£UH  +  6-terms  ,    (77) 


Two  more  simplifications  are  in  order.  First,  as  H  has  a  finite 

jump  across  the  collisional  sublayer,  its  value  is  bounded  inside  the 

1  /? 

sublayer.   Since  the  sublayer  width  is  of  X    order,  the  sublayer 

contribution  to  (77)  can  be  neglected.  We  can  then  use  (42)  for  H  in 
£  >  A  and  (1*8)  for  H  in  £  <  A  without  paying  attention  to  the  sublayer 
solutions.  Secondly,  we  observe  that  the  expressions  in  (42)  and  (48) 
are  independent  of  the  sign  of  u,  and  therefore  we  only  need  to 
calculate  the  integral  in  (79)  for  u  >  0  and  double  the  obtained 
result.  If  we  take  all  these  steps  and  change  the  order  of 
integrations  in  (77),  we  obtain 


'sr 


■S^cVa"1 


v   dv  L. 
1     1 


rA, 


d£<U>   H  + 


rA 

-A 


d£<U>H»]   +   6-terms   ,      (78) 


where  H+  is  given  in  (42)  and  H*  in  (48),  and  L1  is  defined  immediately 
after  (26).  We  then  substitute  (42)  and  (48)  into  (78),  cast  the 
boundary  terms  for  £  «  A»  aside  into  the  6-terms  and  carry  out  the 
integrations.      The   final   result   for    psr    is 


p„_   -   8ir1/26*v  T~1eA1/2q(o)    exp    (-c2/2T)   +    6-terms    , 


sr 


(79) 


where  the  constant  B«  is  given  by 
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e«  -  ^  * 


i 


dn(n-D[n1/2  -  V(n)]/n1/2V(n)    , 
v(n)  -   f      (n  -  cos  e)1/2de/2n   , 


(80) 
(81) 


the  function  q(o)  is  defined  as 


q(o)  -  2  -   [(o+1)/(o+x)]  exp  (-x)dx  , 


(82) 


and  o  -  1/(1+z). 

We  can  now  substitute  (67)  and  (79)  into  (64)  and  set  pg  ■  0  where 
the  6-terms  in  (67)  and  (79)  should  cancel  out.  This  then  gives  us  the 
evolution  equation 


2c2c  3[a  +  YvpA  +  BvpA1/2/e)  ■=  0  , 


or 


A  =  -Yve  A  -  6veA1/2/c  , 


where  6  is 


=  ^TT1/26j((c3/T)q(o)  exp  (-c2/2T)  , 


(83) 


(84) 


(85) 


where  "Y  and  8«  are  positive  constants  given  in   (68)  and  (80), 
respectively,  and  ve  is  also  given  in  (68). 

The  function  q(o)  can  be  expressed  in  terms  of  the  exponential 
integrals  and  has  the  property  3  that  q  >  0  for  z  _<  1 1  and  q  <  0  for 
z  _>  12.  This  shows  that  while  the  isotropic  collision  terms  in  (25) 
involving  the  v-derivat ives  are  stabilizing,  the  pitch-angle  collision 
terms  involving  y-derivatives  are  destabilizing  the  small-v  electrons. 
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Had  we  neglected  the  isotropic  collision  terms,  we  would  have  an 
infinite  nonlinear  growth  rate  and  the  problem  would  become  spurious. 

Using  similar  procedures  we  can  evaluate  the  cosine-component  of 
the  space  charge.  The  lowest  order  result  gives  the  linear  dispersion 
relation,  while  the  higher  order  terms  determine  iji(t)  according  to  eq. 
(8).  We  shall  skip  the  details  of  these  calculations  and  simply  state 
that 


*  -  e6A3/c* 
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VII.   TEMPERATURE  CHANGE  AND  CURRENT  DRIVE 

We  now  calculate  the  change  of  temperature  In  the  background 
Maxwelllan.  In  doing  this,  we  shall  also  be  able  to  demonstrate  the 
cause  of  the  discrepancy  between  the  present  result  and  that  obtained 
in  Ref.   9. 

Let  us  consider  eq.   (23)  and  recall  that  the  e-derivative  of  $  is 

2 
simply  -e  A  sin  6  to  the  lowest  order  in  e.   For  our  purpose  it  is 

convenient  to  rewrite  eq.   (23)  as 


I  f  ■  "Acp=  *  s? 


d6dv(vx-c)F  sin  6  +  yv(Ti~T)  ,  (86) 


where  the  sine  component  ps  of  the  space  charge  is  given  in  (64).  Note 
that  ps  «=  0  according  to  eq.   (7). 

The  integral  in  (86)  can  be  evaluated  using  the  non-resonant 
solution  in  (17)  and  the  resonant  solution  and  the  dominant 
contribution  Is  from  the  A  term  in  (18),  namely 


d6dy(vx-c)F  sin  6  =  2irA/c2  ,  (87) 


which,  upon  substituting  into  (86),  yields 


|  f  -  -e4AA/c2  +  yv(Ti-T)  .  (88) 


h    •      2 
It  can  be  shown  by  using  the  dispersion  relation  that  the  e  AA/c 

o 

term  is  just  the  rate  of  change  of  the  wave  energy  E  /8tt  averaged  over 

6.  Therefore  eq.   (88)  simply  expresses  the  conservation  of  the  sum  of 
the  electron  thermal  energy  and  the  wave  energy  except  for  the  energy 
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transfer  between  electrons  and  ions.   It  also  shows  that  T  is  much 
smaller  than  A  or  v. 

At  this  Juncture,  it  is  of  interest  to  examine  the  result  of 
Zakharov  and  Karpman.  They  used  the  wave  energy  W  ■=  <E  >/8ir  and 
evaluate  its  rate  of  change  according  to 

W  ■=  -  <Ej>  ,  (89) 

where  j  is  the  current  in  the  x  direction. 

In  terms  of  our  non-dimensional  variables,  (89)  can  be  written  as 

e2AA/c2  =  -cAps  -  (A/2tt)  I  d6dv(vx-c)F  sin  6  ,  (90) 

which,  upon  substituting  (83)  and  (87),  becomes 


A  =  -2(A  +  -TveA  +  6veA1/2/e)  +  A  ,  (91  ) 


where  the  last  term  is  due  to  the  integral  in  (90)  and  the  bracketed 
term  comes  from  p  . 

The  procedure  used  in  Ref.  9  is  to  ignore  the  integral  in  (90)  on 
the  ground  that  vx"c  is  small  in  the  nonlinear  layer.  It  further 
neglects  the  first  two  non-resonant  terms  in  the  right  hand  side  of 
(91).  If  we  do  this,  we  would  end  up  with  a  result  that  differs  by  a 
factor  of  two  with  the  B-term  in  (8*0.  The  proper  way  of  determining  A 
is,  of  course,  to  use  ps  -  0  as  in  (83).  If  the  last  term  of  (91)  is 
retained,  Zakharov  and  Karpman  would  obtain  the  correct  result.  In  any 
case,  since  the  last  term  in  (90)  is  equal  to  the  W  term  on  the  left 
hand  side,  the  omission  of  this  term  is  certainly  erroneous. 
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We  further  note  that  the  decay  term  in  Ref.   9  is  of  A-1/2  order 

1  /? 
instead  of  A  '  as  obtained  in  (84).  As  far  as  we  can  determine,  this 

is  related  to  the  first  term  in  the  bracket  of  (75) .   If  this  term  does 

I/O 

not  vanish,  it  would  lead  to  an  A  result.  However  we  have  shown  in 
the  last  section  that  it  must  vanish.  There  are  other  errors  in  Ref. 
9,  such  as  the  sign  of  <f>  in  their  characteristic  variable  e  in  eq. 
(17)  of  Ref.  9  and  the  integration  constant  in  eq.  (2*0  of  Ref.  9. 
These  make  a  complete  comparison  of  the  present  results  with  those  of 
Ref.  9  somewhat  complicated  and  we  shall  not  discuss  the  details 
further. 

Finally  we  wish  to  calculate  the  d.c.  current  j  generated  by  the 
longitudinal  wave 

j  -  NZeve  J  vx(F-FM)d9dv  «  Jn  ♦  jp  ,  (92) 

(93) 
(94) 


Jn  "  NZeve  | 


2tt 

de 

o           J 

n      '    * 

2tt 

de 

0                J 

fr    (F"FM)vxd-v    ' 

Jr  -  NZeve 


where  the  subscripts  n  and  r  denote,  respectively,  the  non-resonant 
region  and  the  resonant  region  contributions  and  the  integration  limits 
are  the  same  as  in  (65)  and  (66). 

The  non-resonant  contribution  j  is  the  usual  quasi-linear  term 
which  can  be  evaluated  by  solving  F,  in  (22)  and  noting  that  only  the 
part  of  F1  that  is  anti-symmetric  in  vx  contributes.  The  resonant 
contribution  Jr  can  be  evaluated  by  substituting  (34),  (39)  and  (42) 
into  (94).   Using  the  symmetry  in  u,  we  obtain  for  jp: 


-38- 

jr  -  NZvee  J  de  J   dv(c+eu)(F-FM) 

-  c2NZveec  [  de  J  dv[c2(H-u2Fc/2T2)  ♦  uG  ♦  u2Fc/t] 

For  large  c,  we  only  need  the  first  two  terms  In  the  bracket  of 
the  last  equation,  and  integrating  by  parts,  we  have 

j_  =   -iJnc3c3NZvpe  I  devdv,    *  dE,  I  u|  (h'-F  /T2) 
r  e   J    1  1 JU  C 

which  can  be  evaluated  by  using  (H2)  and  ( -M 8 )  and  the  final  result  is 


jr  *   1  Ti(e2A/T)3/2NZvee  ,  (95) 


where  6  is  the  constant  given  in  (85).   We  only  note  that  the 
non-resonant  current  j   is  proportional  to  e  while  Jr  -  ej. 
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VIII.   CONCLUDING  REMARKS 

The  main  results  of  the  present  analysis  show  that  the  nonlinear 
decay  term  is  proportional  to  (v/e)c3  exp  (-c  /2T)  while  the 
non-resonant  linear  term  varies  as  v.  This  indicates  that  the 
nonlinear  term  becomes  important  if  the  wave  amplitude  parameter  e  Is 
small.  This  must  be  understood  in  terms  of  the  hypothesis  that  for  the 
present  analysis  to  be  valid  we  must  have  e3  >>  v.  When  the  amplitude 
becomes  really  small  so  that  cJ  <<  v,  the  resonant  layer  is  governed  by 
a  collisional  linear  theory.  In  this  case  the  resonant  layer 
contribution  leads  to  a  linear  Landau  decay  term  that  is  independent  of 
v;  here  v  <<  1  but  v  >>  e3. 

Similarly,  we  note  that  the  resonant  current  is  of  order 
e.Dc->   exp  (-c  /2T),  while  the  non-resonant  current  is  of  e  order. 

Because  of  the  factor  of  c3  exp  (-c  /2T),  the  resonant  terms  will 

1  /  2 
be  important  if  c  is  not  too  large,  for  instance,  c  <  3T   .   In 

1  /  ? 

particular,  if  c  <  21         and  e  is  small,  the  resonant  terms  will 

dominate  and  the  wave  evolution  and  the  current  will  be  primarily 
determined  by  the  resonant  layer  contributions 

A  -  -vA1/2/e  ,     j  -  t3A3/2  .  (92) 

For  the  case  of  c  -  1 ,  the  analysis  can  be  carried  out  with  some 
modifications  which  mainly  involve  the  evaluation  of  the  plasma 
dispersion  function  and  the  treatment  of  the  Fokker-Planck  terms  in  the 
resonant  layer  without  large  c  expansions.  These  modifications  are  not 
expected  to  change  the  qualitative  feature  of  (92). 
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